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Abstract 

We calculate the averaged number v of the lightest supersymmetric particles 
(LSPs) in a shower from the decay of superheavy particle X by generalized DGLAP 
equations. If the primary decayed particles have color charges and the virtuality is 
around 10 13 — 10 14 GeV, the averaged number of the LSPs can become 0(100). As 
the result, the upper limit of the mass of the superheavy particle, whose decay can 
produce the observed abundance of the dark matter, can increase from 10 12 GeV 
to 10 14 GeV. Since the typical scale of the inflaton mass of the chaotic inflation is 
around 10 13 GeV, the decay of the inflaton can produce the observed dark matter 
abundance if the reheating temperature is of order 1 GeV. Even for the standard 
model particles with virtuality Q ~ 10 — 100 TeV, the averaged number of the 
LSPs becomes O(0. 1) for gluon, and 0(0.01) for Higgs, which strongly constrains 
the scenario of non-thermal LSP production from the decay of moduli with 10-100 
TeV mass. 
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1 Introduction 



Recently, the amount of the dark matter (DM) in the universe has been precisely measured 
by WMAP observation [TJ as 

n DM h 2 = 0.1120 ±0.0036. (1.1) 

However, the origin of the DM is still unknown, and it is an important challenge to reveal 
the origin of the DM. One of the most promising candidates for the dark matter is weakly 
interacting massive particle (WIMP), because the thermal production can explain the 
observed amount and such a particle is required to solve the naturalness problem on the 
standard model (SM) Higgs. For example, in the minimal supersymmetric (SUSY) SM 
(MSSM), the lightest SUSY paricle (LSP) becomes stable because of R-parity. And if the 
LSP is a neutralino, it becomes a candidate of the WIMP. We can check this scenario 
if the thermal abundance is calculated by the parameters measured by the experiments. 
However, in principle, it may happen that the calculated thermal abundance is inconsistent 
with the observed value. Moreover, even in the SUSY models, there is a possibility that the 
LSP has no weak interaction like the gravitino or axino. Therefore, it is worth considering 
other possibilities for the production of the LSP. 

In the literature, a lot of studies on the non-thermal production of the LSP have 
been studied [2] -[TO]- The energy density of the LSP x originated from the decay of the 
superheavy field X at the decay of X is estimated as 

Px = m x n x ~ m x vn Xi (1-2) 

where m x , n x , and nx are the mass of the LSP, the number density of the LSP, and the 
number density of the superheavy field X. v is an averaged number of produced LSPs 
from the decay of a single superheavy field X. The parameter v is important to obtain the 
final result, though it is dependent on the explicit models. In many cases, the relation v ~ 
1 is adopted. This is because when the R-parity of the X particle is odd, it is obvious that 
the number of the produced sparticle is around one in the lowest order calculation, and 
when the R-parity is even, the branching ratio to a pair of sparticles becomes the same 
order as the branching ratio to the SM particles in the lowest order calculation due to the 
supersymmetry. However, it has been discussed that the v can be much smaller than 1 
in the special cases. For example, the branching ratio for the decay of the moduli with 
10-100 TeV mass into the gaugino pair can be much smaller than into the gauge boson 
pair because of chirality suppression [4J. In the scenario, the small v ~ 10~ 4 is required to 
obtain the observed abundance of the Wino-like LSP. 

In this paper, we reconsider how many LSPs are produced by one superheavy particle 
decay. Basically, the primary particle in the decay must produce a shower at least if the pri- 
mary particle has color charges. And if the mass of the superheavy particle is much larger 
than the SUSY breaking scale, then it is naturally expected that in the shower, sparticles 
can be produced. We estimate the number of the LSPs in the shower by using super- 
symmetric Dokshitzer-Gribov-Lipatov-Altarelli-Parisi(DGLAP) equations |11[ fT2j [13]. Of 
course, the number of the LSPs in the shower is dependent on the mother MSSM particles 
produced by the decay and on the mass scale of the superheavy particle. We make several 
tables so that our calculation can be applied to various masses of the superheavy particles 
and to various decay modes. 
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In section 2, we briefly remind the scenario on the non-thermal production of the LSP. 
In section 3, we calculate the parameter v numerically. The section 4 is for the discussions 
and summary. 



2 Non thermal production of LSP 

In this section, we estimate the LSP production by the decay of the superheavy particle 
X. We assume that the decay modes are some of the MSSM particles. Supposing that 
the superheavy particle X dominates the energy density of the universe at the decay, the 
radiation energy density p ra d just after the decay can be estimated by the energy density of 
the superheavy field, px = rrix n x, where mx and nx are the mass of X and the number 
density of X, respectively. The reheating temperature T RH can be defined by the relation 
Prad = 9*(Trh)Tr H at the decay. The ratio of the dark matter number density, n x , to the 
entropy density, s, produced by this direct decay process is estimated as 

^direct = u nx. _ v J_fh*A _ 3 h , 9*(T RH ) T rh \\ 

s s m x s 4 g* s {T RH ) m x ^ ' ' 

3, 



= 7.5 x 10--, (jfe) (^) , (2.2) 



where g* and g* s are the number of freedoms of thermalized particles for the radiation 
energy density and for the entropy density, respectively. Usually, the number of produced 
LSPs per one X decay, u, is taken as v ~ 1. In the ratio of the LSP density to the critical 
density, 

^ = ^{^\(^\ms*L\ (2 .3) 



.100 GeV/ V 1MeV / V m x J 
where m x is the mass of the LSP. For the non-thermal production, the reheating temper- 
ature T RH must be larger than 0(1 MeV) not to spoil the success of the Big Bang Nu- 
cleosynthesys (BBN) and must be smaller than the freeze-out temperature of the WIMP 
which is roughly given by the relation Tf reeze _ out ~ m x /20. If the reheating temperature 
is 5 GeV, then the superheavy mass must be taken as 10 12 GeV in order to obtain the 
observed dark matter abundance. 



3 Effect of shower 

If the mass of the superheavy particle is much larger than the SUSY breaking scale, then 
even the showers produced by the primary decay modes include additional SUSY particles. 
If v becomes much larger than the naive expectation, 1, then the above estimation for the 
dark matter abundance must be changed. In this section, we estimate the parameter v by 
using the generalized DGLAP equations [13]. 

Let us briefly outline the physics on this calculation. The primary decay products are 
the MSSM particles and generically off-shell. Since each MSSM particle in the primary 
decay modes has very large virtualities of order mx, it produces a shower. In the shower, 
one virtual particle splits into two other particles with smaller virtualities. When the 
virtuality is larger than the SUSY breaking scale, SUSY particles are also produced in the 
shower. 
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What we would like to know is how many SUSY particles appear in the shower produced 
by the primary MSSM particle / with 0(mx) virtuality. Here, we call the number fi(mx), 
since the number is dependent on the primary MSSM particle /. Then the averaged v is 
written as 

^ = ^Br(X^/)5> 7 (™*), (3-1) 

where Br(X — > f) is the branching ratio of X to the decay mode /. Once we know the 
numbers vi for all / =MSSM particles, we can apply the above relation for the averaged 
v for any models in which the superheavy particle X decays to the MSSM particles. We 
calculated these parameters basically following the technique developed in Refs. [13J. The 
parameters we calculated are listed in Table 1 and 2. 
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Table 1: The number of the LSP, vj, in the parton shower by the MSSM particle / with 
virtuality Q (m x ~ 2Q). We take tan/3 = (H u )/(H d ) = 3.6. 



Let us explain briefly how to calculate these parameters. In order to calculate z/, we 
introduce the fragmentation functions (FFs) Dj(x,Q 2 ) (0 < x < l,m 2 j < Q 2 ), where 
Dj(x,Q 2 ) is the number density of particle J with the energy xQ which are produced in 
the parton shower by the primary field / with initial virtuality Q. Then, the vj is estimated 
by 

Vl (m x )= I dx Dj =LSP (x,m 2 x ). (3.2) 



The fragmentation functions can be obtained as a solution of the DGLAP equation 

^Ln/ M2)= £ E!*m f l d ^ p K Uv)D^/v,Q% (3.3) 



3 



Q[GcV] 


10 11 


10 12 


10 13 


10 14 


10 15 


10 16 




33(34) 


79(86) 


190(220) 


480(590) 


1200(1600) 


3200(4600) 


Ql 


34(35) 


80(86) 


190(220) 


480(590) 


1200(1600) 


3200(4500) 


Qr 


32(33) 


76(83) 


190(220) 


470(580) 


1200(1600) 


3100(4500) 


Qr 


33(34) 


77(84) 


190(220) 


470(580) 


1200(1600) 


01 r\r\ / a f r\ r\\ 

3100(4500) 


k 


1.5(1.6) 


3.4(3.6) 


7.9(9.2) 


r\r\ / r\ \ 

20(26) 


52(78) 


150(250) 


k 


2.5(2.5) 


4.3(4.6) 


8.8(10) 


21(26) 


53(78) 


150(250) 


Ir 


0.20(0.20) 


0.43(0.46) 


1.1(1.2) 


2.8(3.5) 


8.3(11) 


26(40) 


Ir 


1.2(1.2) 


1.4(1.4) 


2.0(2.2) 


3.8(4.5) 


9.3(12) 


27(41) 


9 


66(68) 


160(170) 


390(450) 


980(1200) 


2500(3400) 


6600(9800) 


9 


66(68) 


160(170) 


390(450) 


980(1200) 


2500(3400) 


6600(9600) 


W 


6.6(6.7) 


15(16) 


37(40) 


92(110) 


240(310) 


f~* v r\ / r\c\ t~\\ 

650(920) 


w 


7.3(7.4) 


16(17) 


37(41) 


93(110) 


240(310) 


650(920) 


B 


1.5(1.5) 


3.9(4.0) 


10(11) 


27(30) 


75(88) 


210(260) 


B 


2.4(2.4) 


4.8(4.9) 


11(12) 


28(31) 


k-7 c / r\r\\ 

76(90) 


210(270) 


H u 


2.2(2.2) 


4.7(4.9) 


11(12) 


26(30) 


65(82) 


170(240) 


H u 


3.1(3.1) 


5.6(5.7) 


12(12) 


27(30) 


66(80) 


170(230) 


H d 


1.6(1.6) 


3.5(3.7) 


8.0(9.2) 


20(26) 


52(77) 


150(250) 


H d 


2.5(2.5) 


4.3(4.6) 


8.9(10) 


21(26) 


53(78) 


150(250) 



Table 2: The number of the LSP, vj, in the parton shower by the MSSM particle / with 
virtuality Q {m x ~ 2Q). These are estimated by extraporating the 30(45) data points 
with 10~ 7 < x < 10~ 3 (0.1), which are calculated by generalized DGLAP equations. We 
take tan = 3.6. 



where Pk<-i(x) and oiki{Q 2 ) are splitting functions (SFs)[14j on the three points interac- 
tion which produces K in the MSSM and the running coupling constant for the interaction, 
respectively. As the boundary condition, 

Dj(x, rrij) = Sj8(l — x) (3.4) 

is imposed. This condition is to require that on-shell particle does not produce new particle. 
In order to solve these DGLAP equations, the generalized FFs D^(x, Q 2 ; Ql)(mj < Qo < 
Q) are introduced, which are defined as solutions of DGLAP equations with the boundary 
condition Dj(x,Ql;Ql) = 5j5(l — x). The usual FFs can be written as Dj(x,Q 2 ) = 
Dj(x, Q 2 ; m 2 ). It is quite useful that the usual FFs can be decomposed as 

D J T (x,Q 2 )= f-Df{y,Q 2 -Ql)D J K {xly,Ql), (3.5) 

if Qo is taken as the SUSY breaking scale, the electroweak symmetry breaking scale, or 
hadronization scale, etc. For example, if we take Qo as the SUSY breaking scale msusY, 
then, the virtuality is always larger than the SUSY breaking scale in the DGLAP equations, 
Q > m susYi and therefore, it is sufficient to solve the SUSY DGLAP equation to obtain 
Dj(x, Q 2 \ «t|usy)- I n other words, we can neglect the SUSY breaking effect in calculating 
Dj(x, Q 2 ; t^susy)- Since the above equation can be rewritten as 

/ dxDjix, Q 2 ) = [ d y > ^ ^0) / dzD J K (z, Ql), (3.6) 
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the parameters vj can be obtained as 



vi(m x ) = f dxD J I - LSP (x,m 2 x )= ^ / dy Df{y,m 2 x ;m 2 SVSY ) j dzD 

JO fc'^TMQQA/T J JO 



K (z, ^susy) 



e"gmssm • 



Yl J d v D f(y^ m x\ m lvsY)- ( 3 - 7 ) 

-R"Ssparticles ^ 

Here, in the last equality, we use the assumption that one sparticle with virtuality Q = 
msusY produces only one LSP in the decay, i.e., L dzD J K {z, mg USY ) = 1 f° r an y sparticles 
and the SM particle with virtuality Q = msusY produces no LSP. 

In order to estimate the generalized FFs, we use the program " SHDecay" in which the 
generalized SUSY DGLAP equations are numerically solved[13]. The program "SHDecay" 
includes the all gauge interactions in the SM and the third generation Yukawa couplings. 
We calculated the parameters uj for the third generation fields and the first two generation 
fields separately, but the calculated values are the almost same, so in Table 1 and 2, we 
do not distinguish the third generation fields from the first two generation fields. We 
have several remarks on the calculation. First, the integration in eq. (I3.7p has infra-red 
divergence, because we solved the generalized DGLAP equations in SUSY limit. The 
number of sparticles with smaller energy than their mass becomes quite large. Of course, 
this situation is unphysical. We just introduce the infra-red cutoff for the parameter y as 
y m in = m susy/ m x- (We took msusY = 1 TeV.) Second, in the "SHDecay", the FFs can be 
calculated until y = 10 -7 . This is because of the reliablity of the perturbation. It has been 
noted in Ref.[13] that the energy conservation can be checked in this calculation in the 1 
percent level. Actually, main contribution to the energy comes from the FFs with larger 
y, but for the number of the produced particles the FFs with smaller y is more important. 
For rax < 10 10 GeV, we can calculate the parameter v by integrating the FFs directly 
(see Table 1), but for nix > 10 10 GeV, we have no data. Therefore, we just assume that 
the Dj = Ajy ai , where the aj and Aj are determined by fitting the FFs Dj with 30(45) 
points between 1CT 7 < y < 1CT 3 (0.i), and we can obtain the parameters uj as 

/icr 7 
dyDf(y,m 2 x ;m 2 slJSY ), (3.8) 
i\ fcsparmcies 

where e = msusY/ m x- The results are shown in Table 2. Unfortunately, this approxi- 
mation is not so good especially for large Q, as seen in Table 2 that the estimated values 
for the two different fitting regions have larger discrepancies for larger Q. This is because 
the FFs in smaller x are more important for the estimation of the number of produced 
sparticles in the shower. 

As seen in Table 1 and 2, there is a rough relation between i>j and vj as Vj ~ Vi + 1, 
which is reasonable because the difference of the number of primary sparticle is just one. 

It is interesting that the upper bound of the superheavy mass for obtaining the observed 
DM abundance can be increased from 10 12 GeV to 10 14 GeV, because the inflaton mass for 
the chaotic inflation is around 10 13 GeV. For the superheavy field with the mass around 
10 13 GeV, the ratio of the density can be written as 
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If the reheating temperature of the inflation is around 1 GeV, non-thermal production of 
the LSP can explain the observed DM abundance. 

In the heavy gravitino scenario in which the cosmological moduli problem can be solved 
because the lifetime of the moduli fields becomes shorter than 1 second, the v values for 
Q ~ 10 4 ~ 5 GeV must be important. In Ref. [I], it is noted that in order to produce 
the observed abundance of the Wino-like LSP, v ~ is required. However, in our 

calculation, especially, for the moduli which decays to gluons (Higgs) in certain portion, v 
can be O(0.1)(O(0.01)), which is too large to obtain the observed DM abundance. 



4 Discussion and summary 

In principle, also in the process of the thermalization of the high energy particles, sparticles 
can be produced. Actually, the hard collision between the high energy particle and a 
particle in the thermal bath can produce the sparticles if the center of mass energy is 
larger than the sum of the produced sparticle masses. If the energy of the high energy 
particle, E, is larger than m 2 SUSY /Trh, such processes can be expected. However, as 
discussed in Ref. 0, such production is negligible for the LSP non-thermal production, 
because the thermalization of the high energy particles is so rapid through soft processes 
that the hard collision rarely happens before the thermalization finishes [13]- [1 7] • 

In summary, we calculate the averaged number of the LSPs produced by the decay of 
a single superheavy field X by using generalized DGLAP equations. The number v can be 
0(100) for colored primary particle with the virtuality Q ~ O(10 13 GeV). As the result, 
even by the non-thermal production through the decay of chaotic inflaton with the mass, 
O(10 13 GeV), the observed abundance of the DM can be obtained. Moreover, even if the 
primary decay modes of the heavy moduli with the mass, 0(100 TeV), include only some 
of the standard model particles like gluons or Higgses, the LSPs can produce in the shower, 
which leads to v ~ O(0. 1) or 0(0.01), which is larger than 10 -4 which is required to obtain 
the observed abundance of the LSPs. 

In huge parameter region, the DM produced through the decay of the superheavy field 
X is over produced if the energy density of the superheavy fields dominates the energy 
density of the universe. However, if the produced DM abundance is larger than the thermal 
abundance of the DM, then the pair annihilation process can reduce the abundance. The 
thermal abundance of the DM after the pair annihilation can be estimated as 

n LSP h> ~ 0.25 x (jamjL-Y ( ^ ) (WW) fJL) 4 (4.1) 
VlOOGeW \m 2 LSP (av)J \ T RH ) \g,(T RH ) ) 

by solving the Boltzumann equations [TS]. Since the relation fTL 2 LSP (av) ~ 10 -3 is typical 
for the Wino LSP or Higgsino LSP, then O(100 MeV) reheating temperature can realize 
the observed value for the DM abundance in that case. 

The calculation of the DM abundance in the case in which the energy of the superheavy 
field does not dominate the energy of the universe is straighforward. By multiplying the 

*This value for v ~ 0.1 is reasonable, because there are several processes to produce sparticles through 
virtual gluon, for example, X — > gg* — > gqq, where g and q are gluon and squark, respectively. This 
process is suppressed because of the three body decay, but the number of the final squarks is large, so the 
branching ratio for the process can be large. 
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ratio px/Prad at the decay time to the equations for the DM abundance, we can obtain 
the results. 

We hope our calculation can be applied into many cases in which non-thermal produc- 
tion of DM are taken into account. 
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